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A measurement of the r K tt^Vt branching fraction has been made using 230.2 fb-i of 
data recorded by the BABAR detector at the PEP-II e^e~ collider, located at the Stanford Lin- 
ear Accelerator Center (SLAC), at a center of mass energy ^/s close to 10.58 GeV. We measure 
B{r- ^ K-TY°Ur) = (0.416 ± 0.003 (stat) ± 0.018 (syst))%. 

PACS numbers: 13.35.Dx, 14.60.Fg, ll.30.Hv 



The T is the only lepton with a sufficiently large mass 
to decay to hadrons. Tau decays to hadronic final states 
proceed via W exchange and thus the decay rates to the 
final states containing a strange quark is suppressed by 
the factor (|Kis|/|Kici|)^ relative to the non-strange fi- 
nal states, where \Vud\ and \Vus\ are the moduli of the 
Cabibbo-Kobayashi-Maskawa (CKM) matrix [l], [3| ele- 
ments. For a given value of 3], \Vus\ca.n be deter- 
mined up to unprecendented precision 4, ^ from the 
inclusive sum of the branching fractions of t decays to 
hadronic final states with net strangeness equal to unity. 
This determination of \Vus\ can be perfomed even with- 
out detailed knowledge of the hadronic mass spectrum 
arising due to incompleteness in our understanding of 
some of the intermediate resonance contributions. 

In this paper we present a measurement of the t~ — > 
K^tt^i/t branching fraction^. In recent years, measure- 
ments of the branching fractions for t decays to strange 
hadronic final states have been made using CLEO [6|, 
ALEPH and OPAL [i| data, but have often been lim- 
ited by the size of the available data samples. The high 
luminosity provided by the PEP-II asymmetric-energy 
e+e~ storage rings Q at the Stanford Linear Accelerator 
Center (SLAC) coupled with the large cross-section for 
T+r" -pair production has given us a very large sample 
for studying such decays in the BaBAR detector. 

The BaBAR detector is described in detail else- 
where [lo| . Charged particles are detected and their mo- 
menta measured with a five-layer double-sided silicon ver- 
tex tracker (SVT) and a 40-layer drift chamber (DCH) in- 
side a 1.5T superconducting solenoidal magnet. A ring- 
imaging Cherenkov detector (DIRC) provides additional 
separation power for identification of charged particles 
for momenta greater than 1 GeV/c, and thus comple- 
ments dE/dx measurements in the DCH, useful for the 
identification of charged particles below 1 GeV/c. Ener- 
gies of photons and electrons are measured by a CsI(Tl) 
crystal electromagnetic calorimeter (EMC), and the in- 
strumented magnetic flux return (IFR) is used to identify 
muons. 

The analysis described in this paper is based on a 
data sample corresponding to an integrated luminosity 
of 208.7fb~^ collected at a center-of-mass energy of 
10.58 GeV and 21.5 fb"^ at ^/s = 10.54 GeV. With a 
cross-section for r+r^ pair production averaged over ^/s 



of CT,, = (0.919 ± 0.003) nb [li|, this total data sample of 
230.2 fb"^ contains 211.6 million r+r" pairs. 

Studies of Monte Carlo (MC) simulated events are car- 
ried out for signal and various background samples; r 
pairs are generated with KK2f [l^l and their decays sim- 
ulated with TAUOLA fx^. Signal r decays are modeled 
using form-factors from the Breit-Wigner line-shape of 
K*{892) decays lj|, which nearly saturates the 



K^TT^i/r final state 0, 0, S|- The other r decays into 
any of the possible final states, according to the mea- 
sured branching fractions [l3 |. To estimate non-tau 
backgrounds, samples of T(4S') BB, e+e~ qq 



^ Throughout this paper, the charge conjugate decays are impHed. 



{q = u, d, s, c) and e+e — > fj,'^ fi (7) are generated with 
EvtGen [ll], Jetset7.4 ^ and KK2f [l2| MC programs, 
respectively. The available MC samples are weighted ac- 
cording to their respective size and cross-sections in order 
to match the data integrated luminosity [3] . 

Each event is divided into hemispheres in the center-of- 
mass (CM) frame using the plane perpendicular to the 
thrust axis, which is the direction that maximizes the 
sum of the longitudinal components of the momenta of 
reconstructed particles, both neutral and charged. Only 
events with one charged track in each hemisphere, and 
with both tracks consistent with originating from the 
interaction point (1-1 topology) are selected. The net 
charge of the event is required to be zero. 

To suppress light quark (u, d, s) hadronic events, while 
retaining the relatively large fraction (« 35%) of r lep- 
tons that decay leptonically, we require that one hemi- 
sphere contains a track that is identified either as an elec- 
tron (e-tag) or as a muon (/i-tag). The charged track in 
the opposite hemisphere is then required to be within 
the geometrical acceptance of the DIRC and to be iden- 
tified as a kaon, and inconsistent with an electron. In 
order to reject events where the kaon has decayed or in- 
teracted before reaching the DIRC, a two-dimensional 
requirement on the Cherenkov angle (^c) versus the lab- 
oratory momentum of the candidate kaon (piab) is used: 
dc (rad) < 0.48-1-0.31 *piab (GeV/c). 

Event shape variables are used to discriminate against 
remaining BB and qq backgrounds. The thrust magni- 
tude is required to be greater than 0.9, and the ratio of 
the 2"^^ to the 0**^ Fox- Wolfram moments pL3] is required 
to be greater than 0.5. Also, to account for the substan- 
tial energy carried away by neutrinos in r pair events, 
the total missing momentum in the laboratory frame is 
required to be greater than 0.5 GeV/c. This discriminates 
against Bhabha scattering and //-pair events, as well as 
qq production. Moreover, we remove events in which a 
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decay to two charged pions is identified. 
We further require that the event contains only one 



two separate EMC clusters with an energy above 50 MeV, 
not associated with charged tracks, are considered in the 
analysis. Candidate tt*^ mesons are required to have an 
invariant mass in the range 100 < m^^ < 160 MeV/c^ and 
an energy above 200 MeV in the CM frame. The mass of 
the selected photon pair is constrained to the nominal 7r° 
mass [l^ to improve the purity of the selected tt^s. The 
angle between the charged kaon and tt^ momenta in the 
CM frame is required to be less than 1.0 radian. 

After all the above requirements, there remain 44348 
e-tagged events, and 33764 ^-tagged events. 

The MC simulated events are adjusted to improve their 
accuracy in modeling data events, according to several 
dedicated studies on specific control samples. Charged 
tracks are weighted to compensate for the different par- 
ticle identification (PID) efficiencies between data and 
MC. On average, the MC efficiency is reduced by 15% 
and 3% for muons and kaons, respectively. The electron 
identification is properly simulated and therefore no MC 
efficiency correction is applied. A correction to the tt'^ 
MC efficiency has been obtained from detailed studies 
based on t~ — s- p~i>T and t~ — > TT~i>r events. As a result 
of this study, the tt" MC efficiency is reduced by 2%. 

We estimate the t~ K~h'^Vt selection efficiencies to 
be 1.31%, 0.96% and 2.27% for the e-tag, /.t-tag and com- 
bined samples, respectively, using the signal MC sample 
with all requirements and corrections applied. The effi- 
ciency as a function of the if ~7r° mass is consistent with 
being constant. 

Figure [T] shows the invariant mass spectra of the 
selected K~tt° candidates and simulated backgrounds 
for the combined sample after all the analysis require- 
ments. The contribution from light quarks, cc and 
//+yU~(7) backgrounds is much smaller than the r-pair 
backgrounds. The contribution from BB backgrounds 
is negligible. The K*{%92)'' resonance is seen promi- 
nently above the simulated background. Decays to higher 



K* resonances are expected near 1.4GeV/c^ \\i 
such as the if*{1410)- [IH and K*{\Ai{))- ^,\,ui their 
branching fractions are not well measured yet. These de- 
cays are not included in our simulation of r-pair events, 
but appear to be present in the data near 1.4 GeV/c^. Be- 
low 0.7 GeV/c^, the background is dominated by K~tt'~'tt'~' 
and K'K^tt'-' events, for which the branching fractions 
are only known with relative uncertainties « 40% and 
« 15%, respectively. These uncertainties are taken into 
account in the estimation of the systematic uncertainty 
due to modeling of the r backgrounds. 

The branching fraction B{t~ — > K'^tt'^i't-) is calculated 
as S = [N,,i - iVbkg]/2e,igiV,,, where N^^ = C a^^ = 
2.116 X 10® is the number of produced t+t^ pairs; e^ig is 
the estimated signal selection efficiency; N^^i is the num- 
ber of events selected in data; and TVhv, is the estimated 
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FIG. 1: Distribution of the reconstructed /S'~7r° mass for the 
combined sample in linear (top) and logarithmic (bottom) 
scale. The dots are the data, while the histograms are back- 
ground MC events with selection and efficiency corrections: r 
background (dashed line), qq (dash-dotted line), n'^ fi~ (dot- 
ted line). 



number of background events obtained from MC simula- 
tions. For the combined e-tagged and /i-tagged sample, 
N^^i — 78112 ±280 events and the estimated background 



is A'bkg — 38247 ± 159 events. The branching fraction 
B{t~ K-TT^Vr) is found to be (0.416 ±0.003 (stat))%, 
where the statistical uncertainty comes from the uncer- 



tainty in the number of signal events, N^^ 



Several 



cross checks were performed by determining the branch- 
ing fraction as a function of tag type, 7r° momentum, 
kaon charge and run period; all were found to give con- 
sistent results. 

A summary of the systematic uncertainties is given in 
Table [H The uncertainty in the charged track recon- 
struction efficiency is estimated to be 0.31% per track, 
based on studies of data control samples of r-pair events 
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TABLE I: Relative systematic uncertainties in the r 
K'tt^^Vt analysis. 



Systematic 


e-tag 


M-tag 


Combined 




(%) 


(%) 


(%) 


Tracking efficiency 


0.62 


0.62 


0.62 


■K^ efficiency 


3.26 


3.26 


3.26 


Particle identification 


2.09 


2.34 


2.17 


TT cross section (ctt-t) 


0.31 


0.31 


0.31 


Luminosity [C) 


0.94 


0.94 


0.94 


Signal modeling & MC statistics 


0.38 


0.52 


0.35 


r backgrounds 


L35 


1.35 


1.35 


TOTAL 


4.28 


4.42 


4.32 



decaying to one charged particle on one side and three 
charged particles on the other side. The systematic un- 
certainty associated with the efficiency of detecting a 7r° 
is 3.26%. The stated uncertainty in the charged parti- 
cle identification efficiency represents the combined un- 
certainty for the two charged tracks, (e, K) or (/x, K). 
This uncertainty includes a contribution due to the mis- 
identification of charged pions as kaons. The uncertainty 
associated with the r"'"r~ pair production cross-section is 
0.31% [m and the luminosity determination uncertainty 
is 0.94%. The effects of approximations in the MC signal 
modeling and of the finite MC statistics on the overall 
efficiency are negligible, but have been included in the 
estimation of the systematic uncertainty. The branching 
fractions for several r decay modes that contribute to the 
background, particularly Cabibbo-suppressed decays, are 
not well known. The resulting uncertainty due to the r- 
pair background estimate on the — s- K~7r^Ur branch- 
ing fraction is 1.35%. Backgrounds from other sources 
are very small and their impact of the signal extraction 
is negligible. 

The total systematic uncertainty is the quadratic sum 
of the individual sources described above, and is 4.32% 
for the combined sample. 

In summary, using 211.6 million r"*" r~ pairs recorded 
by the BaBAR detector, we obtain the following result: 

B{t- -> K'TT^I^r) = 

(0.416 ± 0.003 (stat) ± 0.018 (syst)) %. 

This measurement of the branching fraction is the most 
precise to date and is consistent with the existing world 
average, B{t- R-n^Vr) = (0.454 ± 0.030)% 

ACKNOWLEDGEMENTS 

We are grateful for the extraordinary contributions of 
our PEP-II colleagues in achieving the excellent luminos- 
ity and machine conditions that have made this work pos- 
sible. The success of this project also relies critically on 



the expertise and dedication of the computing organiza- 
tions that support BaBAR. The collaborating institutions 
wish to thank SLAC for its support and the kind hospi- 
tality extended to them. This work is supported by the 
US Department of Energy and National Science Foun- 
dation, the Natural Sciences and Engineering Research 
Council (Canada), the Commissariat a I'Energie Atom- 
ique and Institut National de Physique Nucleaire et de 
Physique des Particules (France) , the Bundesministerium 
fiir Bildung und Forschung and Deutsche Forschungsge- 
meinschaft (Germany), the Istituto Nazionale di Fisica 
Nucleare (Italy), the Foundation for Fundamental Re- 
search on Matter (The Netherlands), the Research Coun- 
cil of Norway, the Ministry of Science and Technology of 
the Russian Federation, Ministerio de Educacion y Cien- 
cia (Spain), and the Science and Technology Facilities 
Council (United Kingdom). Individuals have received 
support from the Marie-Curie IFF program (European 
Union) and the A. P. Sloan Foundation. 



* Deceased 

t Now at Tel Aviv University, Tel Aviv, 69978, Israel 

* Also with Universita di Perugia, Dipartimento di Fisica, 
Perugia, Italy 

^ Also with Universita della Basilicata, Potenza, Italy 
^ Also with Universitat de Barcelona, Facultat de Fisica, 
Departament ECM, E-08028 Barcelona, Spain 
[1] N. Cabibbo, Phys. Rev. Lett. 10, 531 (1963). 
[2] M. Kobayashi and T. Maskawa, Prog. Theor. Phys. 49, 
652 (1973). 

[3] M. Jamin, J. A. Oiler, and A. Pich, Phys. Rev. D 74, 
074009 (2006). 

[4] E. Gamiz, M. Jamin, A. Pich, J. Prades, and F. Schwab, 

JHEP 01, 060 (2003). 
[5] K. Maltman and C. E. Wolfe (2007), arXiv:hep- 

ph/0703314. 

[6] M. Battle et al. (CLEO), Phys. Rev. Lett. 73, 1079 
(1994). 

[7] R. Barate et al. (ALEPH), Eur. Phys. J. C 10, 1 (1999). 
[8] G. Abbiendi et al. (OPAL), Eur. Phys. J. C 35, 437 
(2004). 

[9] PEP-II collaboration (1993), SLAC-418. 
[10] B. Aubert et al. (BABAR), Nucl. lustrum. Meth. A 479, 
1 (2002). 

[11] S. Banerjee, B. Pietrzyk, J. M. Roney, and Z. Was, 

arXiv:0706.3235 [hep-ph]. 
[12] B. F. L. Ward, S. Jadach, and Z. Was, Nucl. Phys. Proc. 

Suppl. 116, 73 (2003). 
[13] S. Jadach, Z. Was, R. Decker, and J. H. Kiihn, Comput. 

Phys. Commun. 76, 361 (1993). 
[14] W.-M. Yao. et al. (Particle Data Group), Journal of 

Physics G 33, 1 (2006). 
[15] D. J. Lange, Nucl. Instrum. Meth. A 462, 152 (2001). 
[16] T. Sjostrand, Comput. Phys. Commun. 82, 74 (1994). 
[17] G. C. Fox and S. Wolfram, Nucl. Phys. B149, 413 (1979). 
[18] M. Finkemeier and E. Mirkes, Z. Phys. C 72, 619 (1996). 
[19] K. Maltman, Nucl. Phys. Proc. Suppl. 123, 149 (2003). 
[20] E. Gamiz, M. Jamin, A. Pich, J. Prades, and F. Schwab, 



7 



Phys. Rev. Lett. 94, 011803 (2005). (1999). 
[21] R. Barate et al. (ALEPH), Eur. Phys. J. C 11, 599 



